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INTRODUCTION 



One of the main problems that optical quantum computing has to overcome is the efficient construction of 
two-photon gates [1] . We can use Kerr nonlinearities to induce a phase shift in one mode that depends on the 
photon number in the other mode, and this nonlinearity is sufficient to generate a universal set of gates [2]. 
However, passive Kerr media have typically small nonlinearities (of the order of 1 CP 16 cm 2 sV^ 1 [3]). We can 
also construct large Kerr nonlinearities using slow light, but these techniques are experimentally difficult [4]. 
On the other hand, we can employ linear optics with projective measurements. The benefit is that linear 
iy) ' optical schemes are experimentally much easier to implement than Kerr-media approaches, but the downside is 
that the measurement-induced nonlinearities arc less versatile and the success rate can be quite low (especially 
CO ' when inefficient detectors are involved). However, Knill, Laflamme and Milburn [5] showed that with sufficient 
ancilla systems, these linear-optical quantum computing (LOQC) devices can be made near-deterministic with 
■ only polynomial resources. This makes linear optics a viable candidate for quantum computing. Indeed, many 
C"> { linear optical schemes and approaches have been proposed since [6-11], and significant experimental progress 
, l t . has already been made [12, 13]. 

The general working of a device that implements linear optical processing with projective measurements is 
shown in Fig. 1. The computational input and the ancilla systems add up to N optical modes that are subjected 
to a unitary transformation U , which is implemented with beam splitters, phase shifters, etc. This is called 
an optical iV-port device. In order to induce a transformation of interest on the computational input, the 
output is conditioned on a particular measurement outcome of the ancilla system. For example, one can build 
a single-photon quantum nondemolition detector with an optical iV-port device [14]. In general, iV-port devices 
have been studied in a variety of applications [15]. 

The class of such devices of interest here is that in which a unitary evolution on the computational input 
is effected. To date these devices have been proposed and studied on a more-or-less case by case basis. Our 
approach is to address this class in a more general way, and identify the conditions that such a device must 
satisfy to implement a unitary evolution on the computational input. Once that unitary evolution is established, 
an effective photon nonlinearity associated with the device can be identified. 

In this paper, we present necessary and sufficient conditions for the unitarity of the optical transformation 
of the computational input, and we derive the effective nonlinearities that are associated with some of the 
more common optical gates in LOQC. We begin section II by introducing the formalism. In sections III-V, 
we examine the transformation equation under the assumption that it is unitary. We show that there are two 
necessary and sufficient conditions for the transformation to be unitary and we provide a simple test condition. 
In section VI, we expand the formalism and conditions to include measurement dependent output processing 
(see Fig. 2), which is used in several schemes. In section VII, we show how the formalism can be applied 
to quantum computing gates. We choose as examples two quantum gates already proposed, the conditional 
sign flip of Knill, Laflamme, and Milburn [5], and the polarization-encoded CNOT of Pittman et al. [7]. Our 
concluding remarks are presented in section VIII, where we note that our main results extend to devices where 
the unitary transformation U is more general than those implementable with linear optics alone. 
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II. THE GENERAL FORMALISM 



We consider a class of optical devices that map the computational input state onto an output state, conditioned 
on a particular measurement outcome of an ancilla state (see Fig. 1). We introduce a factorization of the entire 
Hilbert space into a space He involving the input computing channels (i.e., both "target" and "control" in a 
typical quantum gate), and a Hilbert space Ha involving the input ancilla channels, 

H = He ® U A • 

We assume that the input computing and ancilla channels arc uncorrelated and uncntangled, so we can write 
the full initial density operator as p® a, where p is the initial density operator for the computing channels, and 
a the initial density operator for the ancilla channels. 

Let U be the unitary operator describing the pre-measurement evolution of the optical multi-port device. At 
the end of this process we have a full density operator given by U (p® a) UK In anticipation of the projective 
measurement, it is useful to introduce a new factorization of the full Hilbert space into an output computing 
space H c and a new ancilla space Ha, 

H = He <8> Ha ■ 

The Von Neumann projective measurements of interest are described by projector- valued measures (or PVMs) 
of the type {P, i — P}, where I is the identity operator for the whole Hilbert space, and the projector P is of 
the form 

P = I c ^s- k \k)(k\ , (1) 

k 

where I c is the identity operator in Hq, and we use Roman letters with an overbar, e.g., Ifc), to label a set of 
orthonormal states, (k\l) — spanning the Hilbert space Ha', each is equal to zero or unity. The number 
of nonzero sj, identifies the rank of the projector P in Ha- "Success" is defined as a measurement outcome 
associated with the projector P, and the probability of success is thus 

d(j>) = Trc,A(U(p<2>a)Ulp) . (2) 

Clearly, in general d(p) depends on the ancilla density operator a, the unitary evolution U , and the projector 
P, as well as on p. However, we consider the first three of these quantities fixed by the protocol of interest 
and thus only display the dependence of the success probability on the input density operator p. In the event 
of a successful measurement, the output of the channels associated with Hq is identified as the computational 
result, and it is described by the reduced density operator 

__ Tr^ (PU(p®<j)WP) 
P Tr CA (U{p®a)WPy 

For any p with d(p) ^ 0, this defines a so-called completely positive (CP), trace preserving map T that takes 
each p to its associated p: p = T(p), relating density operators in He to density operators in Hq- It will be 
convenient to write T{p) = V(p)/d(p), where 

V(p) = Tr A (PU (p <g> a) C/ 1 P) (4) 

is a linear (non-trace preserving) CP map of density operators in He to positive operators in H c that is defined 
for all density operators p in He- We restrict ourselves to density operators p over a subspace Sc of He- This 
is usually the subspace in which the quantum gate operates. 

As an example, consider the gate that turns the computational basis into the Bell basis. In terms of polar- 
ization states, the subspace Sc might be spanned by the computational basis {\H, H), \H, V), \ V, H), \V, V)} 
(whereas He is spanned by the full Fock basis). The Bell basis on Sc is then given by {|* + ), I* - ), |$ + ), l^ - )}, 
where 

= -L (\H, V) ± \V,H)) and |d>±) = -L (\H, H) ± \V, V)) . 

This gate is very important in quantum information theory, because it produces maximal entanglement, and its 
inverse can be used to perform Bell measurements. Both functions are necessary in, e.g., quantum teleportation 
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[16]. However, it is well known that such gates cannot be constructed deterministically, and we therefore need 
to include an ancilla state a and a projective measurement. We consider gates such as these in this paper. 
Suppose the subspace Sc is spanned by a set of vectors labeled by Greek letters, e.g., \a). We can then write 

p = ^|a)p^(/3|, (5) 

where p af3 = (a\p\(3). We identify a convex decomposition of the ancilla density operator a as 

o- = ^Pi \Xi) (Xi\ , 

i 

where the normalized (but not necessarily orthogonal) vectors are elements of Ha, and the pi are all 
non-negative and sum to unity, 

i 

We can then use (3) to write down an expression for the matrix elements of p. Note that it is possible to 
work with the eigenkets of a so that {|Xi)} is an orthonormal set; however, this does not simplify the analysis 
so we do not introduce the restriction. Furthermore, dealing with non-orthogonal states in the ancilla convex 
decomposition may be more convenient, depending on the system of interest. Choosing an orthonormal basis 
of H c that we label by Greek letters with overbars, e.g., \a), we find 

^ fi, = EE( w ffo»))^( w So»))*. ( 6 ) 

0,7 i,k 

where 

w ffo») = -*^«*l< a i) tr (^i»»- 

Note that 

a i,k 

which is confirmed by 

Tr C (p) = E ~P && =Y.P" P = Tr cW. (7) 

a 

This last equation follows immediately from (3), since T is a trace preserving CP map and Trc(p) = 1. 

In this paper, we consider a special class of maps that constitute a unitary transformation on the computa- 
tional subspace Sc- In particular, such transformations include the CNOT, the C-SIGN, and the controlled bit 
flip. These are not the only useful maps in linear optical quantum computing, but they arguably constitute the 
most important class. Before we continue, we introduce the following definition: 

Definition: We call a CP map p —* p = T{p) an operationally unitary transformation on density operators p 
over a subspace Sc if and only if: 

• For each p over the subspace Sc we have d(p) ^ 0, and 

• For each p defined by Eq. (5) over the subspace Sc, the map T(p) yields a p given by 

p = Y J \^)p aP {M, (8) 

a,0 

where the \v a ) are fixed vectors in Tic satisfying {u a \9p) = (a\[3) = 5 a p. 

This forms the obvious generalization of usual unitary evolution, since it maintains the inner products of 
vectors under the transformation. Much of our concern in this paper is in identifying the necessary and sufficient 
conditions for a general map T{p) of Eqs. (3) and (6) to constitute an operationally unitary map. We begin in 
the next section by considering what can be said about such maps. 
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III. CONSEQUENCES OF OPERATIONAL UNITARITY 



In this section we restrict ourselves to CP maps T(p) that are operationally unitary [see Eqs. (5) and (8)] for 
density operators p over a subspace Sc of He- The linearity of such maps implies that the convex sum of two 
density operators is again a density operator: 

pc = xp a + (1 - x)p b , 

with < x < 1. Applying Eqs. (5) and (8) to the three density operators p a , Pb, and p c it follows immediately 
that 

Pc = Xp a + (1 - x)p h . (9) 

Now a second expression for p c can be worked out by using the defining relation (3) directly, 

* " ™ " ^7 (10) 

xV{p a ) + (1 - x)V{ Pb ) 
xd(p a ) + (1 - x)d(pb) 
xd(p a )p a + (1 - x)d(p b )pb 
xd(p a ) + (1 - x)d(pb) 

where in the second line we have used the linearity of V(p) (4) and d(p) (2), and in the third line we have used 
the corresponding relations for p a in terms of p a , and pb in terms of p b . Setting the right-hand-sides of Eqs. (9) 
and (10) equal, we find 

x(l - x) [d(p b ) - d(p a )} (p a - p b ) = 0. (11) 

Since it is easy to see from Eqs. (5) and (8) that if p a and pb are distinct then p a and pb are as well; choosing 
< x < 1 it is clear that the only way the operator equation (11) can be satisfied is if d(p a ) = d(pb). But since 
this must hold for any two density operators acting over Sc, we have established that: 

• If a map T(p) is operationally unitary for p (acting on a subspace Sc), then d{p) is independent of p: 
d(p) = d, for all p acting on that subspace. 

With this result in hand we can simplify Eq. (6) for a map that is operationally unitary, writing 

^ = EE^V 7 B 7 )*> ("J 

0,1 j 

where now 

wf = w^.=s k] l^((k\(a\)U(\(3)\xi)) 

is independent of p; we have also introduced a single label J to refer to the pair of indices k,i. A further 
simplification arises because the condition of operational unitarity guarantees that the subspace S c of H c , over 
which the range of density operators p generated by T{p) act as p ranges over Sc, has the same dimension as 
Sc- We can thus adopt a set of orthonormal vectors \a) that span that subspace S c , and the matrices w"' 5 are 
square. 

At this point we can formally construct a unitary map on Sc- p = U(p), which is isomorphic in its effect on 
density operators p with our operationally unitary map p = T(p). We do this by associating each \a) with the 
corresponding \a), introducing a density operator p acting over Sc, and putting 

~ct5 —OiS 



p ad = p ad , (13) 
Mf = wf. 



The unitary map p = U(p) is defined by the CP map 



0,1 j 
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or simply 

p = Y,Mj P M\. (14) 
j 

This is often what is done implicitly when describing an operationally unitary map, and we will see examples 
later in section VII; here we find this strategy useful to simplify our reasoning below. 
Since the map p = U{p) is unitary it can be implemented by a unitary operator M, 

p = MpM\ 

where Aft = M' 1 . Thus (Mi, M 2 , ....) and (M, 0, 0, ....), where we add enough copies of the zero operator so 
that the two lists have the same number of elements, constitute two sets of Kraus operators that implement the 
same map p = U{p). From Nielsen and Chuang [17] we have the following theorem: 

Theorem: Suppose {Ei, . . . , E n } and {Fi, . . . , F m } are Kraus operators giving rise to CP linear maps £ and 
T respectively. By appending zero operators to the shorter list of elements we may ensure that to = n. 
Then £ = T if and only if there exists complex numbers Ujk such that Ej — VJ fe u jkPk, and Uj k is an 
to x to unitary matrix. 

Hence, (Mi,M2,....) must be related to (M, 0, 0,....) by a unitary matrix, and each Mj is proportional to 
the single operator M. This proof carries over immediately to the operationally unitary map T(p) under 
consideration, and we have 

• If a map T(p) is operationally unitary for p acting over a subspace Sc, then for fixed k and i the square 
matrix defined by 

w Z = s k^ ((k\(a\)um\ Xi )), 

with a labeling the row and (i the column, either vanishes or is proportional to all other nonvanishing ma- 
trices identified by different k and i. We can thus define a matrix w a/3 proportional to all the nonvanishing 
w^. such that we can write our map (12) as 

p^ = ^«, a V 7 (^ 7 )*- (15) 

It is in fact easy to show that the two necessary conditions we have established here for a map T{p) to be 
an operationally unitary transformation are also sufficient conditions to guarantee that it is. We show this in 
section V. First, however, we establish a simple way of identifying whether or not d(p) is independent of p. 

IV. THE TEST CONDITION 

In this section we consider a general map T(p) of the form of Eq. (3), and seek a simple condition equivalent 
to the independence of d(p) on p for all p acting over Sc- To do this we write d(p) of Eq. (2) by taking the 
complete trace over He and Ha rather than over H c and Hj, 

d(p) - Tr c ,A(U(p®tT)U*P) 
= Tr c ,A((p®<T)U*PU) 
= Tr c (pT) 

where we have introduced a test operator T over the Hilbert space He as 

T = Tt a (a U^PU) , 

which does not depend on p. The operator T is clearly Hermitian; it is also a positive operator, since the 
probability for success d(p) > for all p. We can now identify a condition for d(p) to be independent of p: 
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Theorem: d(p) is independent of p, for density operators p acting over a subspace Sc of He, if and only if 
the test operator T is proportional to the identity operator Is c over the subspace Sc- We refer to this 
condition on T as the test condition. 

Proof: The sufficiency of the test condition for a d{p) independent of p is clear. Necessity is easily established 
by contradiction: Suppose that d{p) were independent of p but T not proportional to I$ c ■ Then at least 
two of the eigenkets of T must have different eigenvalues; call those eigenkets \p a ) and \pb)- It follows 
that d(p a ) 7^ d(pb), where p a = \p a ) (p a \ and pb = \pb) (Hb\, in contradiction with our assumption. □ 

When the test condition is satisfied we denote the single eigenvalue of T over Sc as r, i.e., T = tIs c - Then 
d{p) — t, and r is identified as the probability that the measurement indicated success. For any given protocol 
the calculation of the operator T gives an easy way to identify whether or not d(p) is independent of p. 

V. NECESSARY AND SUFFICIENT CONDITIONS 

We can now identify necessary and sufficient conditions for a map p = T(p), to be an operationally unitary 
map for p acting on a subspace Sc of Tic- They are: 

1. The test condition is satisfied: Namely, the operator 

T = Ty a (a U f PU) 
is proportional to the identity operator Is c over the subspace Sc- 

2. Each matrix 

™Z = s- k ^((k\(a\)U(\/3)\ Xi )), 

identified by the indices k and i, with row and column labels a and [3 respectively, either vanishes or is 
proportional to all other such nonvanishing matrices; here r is the eigenvalue of T . 

The necessity of the first condition follows because it is equivalent to the independence of d(p) on p, which 
was established above as a necessary condition for the transformation to be operationally unitary, as was the 
second condition given here. So we need only demonstrate sufficiency, which follows immediately: If the first 
condition is satisfied then d(p) — r is independent of p, and if the second is satisfied then, from Eq. (12), we 
can introduce a single matrix w a/3 such that (15) is satisfied. Then 

]T p && = ]T t^i ]T {w a -')* w & v. 

a /3,7 a 

Now the Hcrmitian matrix 

a 

must in fact be the unit matrix: F 7/9 = 5 7) g, otherwise we would not have 

a (3 

for an arbitrary p over Sc, and we know our general map p = T(p) satisfies that condition [see Eq. (7)]. Thus 
w a P is a unitary matrix, and from the form of Eq. (15) of the map from p to p it follows immediately that the 
map is operationally unitary [see Eqs. (5) and (8)]. 

The physics of the two necessary and sufficient conditions given above is intuitively clear, and indeed the 
results we have derived here could have been guessed beforehand. For if the probability for success d(p) of the 
measurement were dependent of the input density operator p, by monitoring the success rate in an assembly 
of experiments all characterized by the same input p, one could learn something about p, and we would not 
expect operationally unitary evolution in the presence of this kind of gain of information. And the independence 
of the nonvanishing matrices w?? on k and i, except for overall factors, can be understood as preventing the 
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'mixedness' of both the input ancilla state a and the generally high rank projector P, from degrading the 
operationally unitary transformation and leading to a decrease in purity. 

If a map is found to be operationally unitary, we can introduce the formally equivalent unitary operator M 
on He, as in Eq. (13), which can then be written in terms of an effective action operator Q, 

M = e~ iQ / h . (16) 

The operator Q can be determined simply by diagonalizing M, and its form reveals the nature of the Hamiltonian 
evolution simulated by the conditional measurement process. We can define an effective Hamiltonian -ff e // that 
characterizes an effective photon nonlincarity acting through a time t e ff by putting H e ff = Q/t e ff, where t e ff 
can be taken as the time of operation of the device. 

In a special but common case, the input ancilla state is pure and the projector P is of unit rank in H A . For 
cases such as this there is only one matrix w a @ in the problem, and thus there is only a single necessary and 
sufficient condition for the map to be operationally unitary: 

• In the special case of a projector P of rank 1 in H&, where P — I c ® |-K") (K\, and a pure input ancilla 
state, a = \x) (xl> then map p = T{p) is operationally unitary for p acting on a subspace Sc of He if and 
only if T satisfies the test condition. Here 

T=(x\U^PU\x), 

which is an operator in He- If it does satisfy this condition, then the transformation is given by 

p &s = w &f) p f) ~< (w E ^ * , (17) 
0,1 

where 

™ &P = \fl((K\ (<*\)U(\® 
and r is the single eigenvalue of T over Sc- 

VI. GENERALIZATION TO INCLUDE FEED-FORWARD PROCESSING 

Suppose that the measurement outcome of the ancilla does not yield the desired result, but that it signals that 
the output can be transformed by simply applying a (deterministic) unitary mode transformation on the output 
(see Fig. 2). This is called feed- forward processing and is widely used. For example, in teleportation, Alice 
sends Bob a classical message which allows him to correct for 'wrong' outcomes of Alice's Bell measurement. 
Here, we can explicitly take into account feed-forward processing. 

Suppose the projective measurement is characterized by a set of projectors, each identifying a different 
detection signature, {Pm, P(2), ■■■P(n), -Pi}; where 

N 

p± = i-Y, p w> 

L=l 

and 

P(L)=I C ®^s Lrk \k)(k\. 
k 

All the s L £ are equal to zero or unity, such that 

P(L)P(L>) = P(L)$LL>- 

Here success arises if the measurement outcome is associated with any of the operators P(l)- And if outcome L 
is achieved, then the computational output is processed by application of the unitary operator Vj^) acting over 
H c . The probability of achieving outcome L is 



d (L) (p) = Tr CtA (U(p®a)U^P (L) ) 
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and if outcome L is achieved the feed-forward processed computational output is then 

_ %) [Yr A (P (L) U{p®cj) WP (L) )] V{ L) 

P{L) ~ Tr c ,A(U(p®v)WP (L) ) 

which defines a map pi L \ = Til){p) f° r those p for which diL){p) 0- I n this more general case we define the 
set of maps \T(l)} to be operationally unitary for density operators p over the subspace Sc when: 

• For each p over the subspace Sc at least one of the d<L){p) =h 0, and 

• For each p over the subspace Sc, for each L for which d^(p) ^ the map T(l)(p) yields a p^ of the 
form of Eq. (8), independent of L. 

The kind of arguments we have presented above can be extended to show that the necessary and sufficient 
conditions for such a set of maps to be operationally unitary for density operators p over the subspace Sc are: 

1. Test conditions are satisfied: The operators 

T (L) = Ti A (aU*P {L) U) 

are each proportional to the identity operator Isc over the subspace Sc- The proportionality constants 
T(£) need not be the same for all L. 

2. Omitting matrices associated with any L for which T( £ ) = 0, each matrix 

<m = s ^\[t^^ v $ «*l ^ u m ' 

identified by the indices L, k, and i, with row and column labels a and f3 respectively, either vanishes or 
is proportional to all other such nonvanishing matrices. 

The probability of success is ^ L TVn = t. This expanded formalism applies to the feed- forward schemes 
discussed by Pittman et al. [13] and the teleportation schemes of Gottesman and Chuang [6]. In devices such 
as these, a measurement provides classical information that is used in the subsequent evolution of the output 
state. 

In a common special case, the input ancilla state is pure, a = \x) (xl> an d each of the projectors P(l) is of 
unit rank in H A , P^ = l c ® \ki) Here the two necessary and sufficient conditions for the set of maps to 

be operationally unitary for density operators p over the subspace Sc simplify to: 

1. All the operators 

T {L ) = (x\U*P { l)U\x) 

over He satisfy the test condition. 

2. Omitting matrices associated with any L for which T( L ) = 0, each matrix 

V (L) - x 

identified by the indices L, with row and column labels a and (3 respectively, either vanishes or is propor- 
tional to all other such nonvanishing matrices. 

If these conditions are met, then the operationally unitary transformation is given by 

0,1 

which is independent of L. 

Another extension of the standard Von Neumann, or projection, measurements is to the class of measurements 
described by more general positive operator valued measures, or POVMs. These can be used to describe 
more complicated measurements, often resulting from imperfections in a designed PVM. Our analysis can be 
generalized to POVMs by expanding the ancilla space, and then describing the POVMs by PVMs in this 
expanded space. In some instances operationally unitarity might still be possible; in others, the extension 
would allow us to study of the effect of realistic limitations such as detector loss and the lack of single-photon 
resolution. 
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VII. EXAMPLES 



In this section we will apply the formalism developed above to two proposed optical quantum gates for 
LOQC. The straightforward calculation of the effects of these gates presented in the original publications make 
it clear that they are operationally unitary; our purpose here is merely to illustrate how the approach we have 
introduced here is applied. 

To evaluate the test operators and matrix elements w^ 3 ^ it is useful to have expression for quantities 

such as UciqU^ , where we use capital Greek letters as subscripts on the letter a to denote annihilation operators 
for input (computing and ancilla) channels; similarly, we use a A to denote annihilation operators for output 
(computing and ancilla) channels. We now characterize the unitary transformation U by a set of quantities U 
* 2A that give the complex amplitude for an output photon in mode A given an input photon in mode ft. That 
is, 

U (4 |vac)) = £ U* A (4 |vac>) , (18) 

A 

where |vac) is the vacuum of the full Hilbert space H. Since only linear optical elements are involved we have 
f/t | vac) = | vac), and it further follows from (18) that 

Ua\^ = % A a A' ( 19 ) 

A 

or 

Ua n tf = Y^U nA a A . (20) 
A 

Using the commutation relations satisfied by the creation and annihilation operators, it immediately follows that 
the matrix Uq A , which identifies the unitary transformation U, is itself a unitary matrix. Certain calculations 
can be simplified by its diagonalization, but for the kind of analysis of few photon states that we require this is 
not necessary. We will need to express, in terms of few photon states with respect to the decomposition Hq ® 
Ha, the result of acting with U on few photon states of the decomposition He ®Ha', this follows directly from 
(19). For example, denoting by |lni2n 2 ) the state with one photon in mode Oi and two in mode O2, we have 

tf|l ni 2n 2 ) = ^k^Hvac) (21) 
= ^ (U^rf) (Ualtf) (t/a^f/t) |vac) 

= 4 E Uq^U^a^As ( a L a A 2 a A 3 ) l vac >' 

V Ai,A 2 ,A 3 

and doing the sums in the last line allow us to indeed accomplish our goal. 



A. KLM conditional sign flip 



The first example we consider is the conditional sign flip discussed by Knill, Laflamme, and Milburn [5]. Note 
that in this case the input ancilla state is pure, there is no feed-forward processing, and the projector P is of 
unit rank in Hx- The necessary and sufficient conditions for the transformation to be operationally unitary are 
those of the special case discussed in section V. The gate consists of one computational input port (labeled 1) 
and two ancilla input ports (2 and 3). The projective measurement is performed on two output ports (b,c) and 
the one remaining port is the computational output (a). The subspace So is spanned by the Fock states |0), 
|1), and 1 2) in each optical mode. 

The pre-measurement evolution, which is done via beam splitters and a phase shifter, is given by the unitary 
transformation U and characterized by the matrix 
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U = U" 



The ancilla input state is 



1-V2 2" 1 / 4 (3/V2-2) 1 / 2 

1/2 1/2-1/V^ 
(3/V2-2) 1 / 2 1/2-1/V2 V2-1/2 



Ix) = 4 l vac A) 



(22) 



(23) 



denoting a single photon in the 2 mode, where |vacA) denotes the vacuum of Ha- The projective measurement 
operator is given by 

P = J c ® |.K") = Iq ® a\ |vac^) (vac^| a,b 

which corresponds to the detection of one and only one photon in mode b, and zero photons in mode c. The 
basis states that define the subspace Sc are 



|0) = |vac c ), |1) =4|vac c ), |2) 



and the basis states of Hq are 



|0) = |vac c ), 1 1) = 4 |vac c ) , |2) 



(A 

V2 



V2 



|vac c ) 



|vac c ) 



In order to evaluate the test function, we first write 

U^PU = J2 U1 ( a l \™ c a) ® l«>) «a| ® (vac^l a h ) U 

a 

and look at the matrix elements 



(M®<x|)tftPtf(| X >®|/J>) (24) 
= ^2 ® ( VSlC a\ 02) (a\ |vac^) ® |a)^ ((a| ® (vac^| a b ) {/ (4 l vac /i) ® 

over the computational subspace, 5c- The calculation is straightforward. Applying the operator U on each of 
the states a\ Ivac^) ® \(3) gives the following states in the Hq ® 7^ decomposition 



c/ (4 l vac ^) ® 

U (a\ \vac A } ® |1)) 



(2- 1/4 4 + ^4 + 



2- 1/4 4 + \a\ + 



1 1 

2 ~ 71. 

1 1 

2 ~ 71. 



oT ) |vac) 



1 - V21 4 + 2- 1 / 4 4 



V2 



1/2 



4 l vac > 



c/ (4 ivac.4) ® 12)) = -L (Vv 4 4 + 14 + [I _ _L] 4) 



1 - y/2 4 + 2-^ 4 4 + 



V2 



- 2 



1/2 



I vac) 
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and we can then separately evaluate the terms in the sum (24), noting that the non-zero elements are 

2 1 



((0| <g> (vac A \a b ) U (a\ \va,c A ) ® |0)) 
((l| ® (vac^| db) U (a\ \v&ca) ® 
«2~| <g> (vac^| a b ) U (a\ \v&c A ) ® |2)) 



The test operator T is then 



( a i) 2 (aif 
|vac c ) (vac c | + a} |vac c ) (vac c | Oi + |vac c ) (vac c | — ^=r- 



= ^sc, 

and is indeed a multiple of the unit operator in the computational input space. The probability of a success- 
indicating measurement is 1/4, independent of the computational input state. Since this test condition is 
satisfied, the transformation (17) is operationally unitary. The terms of the transformation matrix w al3 can be 
calculated noting that the non-zero {R\ (a\ U \0) \\) terms are 



(ir|(o|t/|o)| x ) = 

(K\(T\U\l)\ X ) = 



(K\(2\U\2)\ X ) = 



1 



and since r = 1/4 the non-zero elements of the transformation matrix are 

w°° - 1, 
w Tl = 1. 



w 



22 



which corresponds to the conditional sign flip, since with probability 1/4 the gate takes the input state \ip) = 
ceo |0) + ax |1) + oi2 |2) and produces the state |i/>) = cuo |0) + ol\ |l) — tti |2). 

This map can be seen to exhibit an effective nonlinear interaction between the photons, since the formally 
equivalent unitary map (see section III) is characterized by the unitary operator M (13), 

i>) = oto |0) + ay |1) - a 2 |2) = M (a |0) + ai |1) + a 2 |2)) , 

which can be written in terms of an effective action operator Q (16), where we can take 



Q 



'). 



with n the photon number operator. But such an effective action operator exists only if we restrict ourselves to 
the three-dimensional subspace Sc, spanned by the kets |0), |1), and |2). For consider an attempt to expand this 
subspacc to that spanned by the kets (|0) , |1) , |2) , |3)). The device guarantees that a computational input of 
three photons can only produce a computational thrcc-photon output, since a successful measurement requires 
the detection of one and only one photon in the ancilla space. The test operator is therefore still diagonal in 
the photon number basis. However, we find 



((3| (E> (vac^j a b ) U jvac^) ® |3)^ 



and thus the test operator T is no longer a multiple of the unit operator in this enlarged subspace. In this 
larger space the probability of a success-indicating measurement is dependent on the input, and the map is not 
operationally unitary. 



12 



B. Polarization encoded CNOT 

The second example is the polarization-encoded Gottesman-Chuang protocol discussed by Pittman et al. [7]. 
In this case the input ancilla state is pure, there is feed-forward processing, and there are several projectors 
P(i) of unit rank in H^. The necessary and sufficient conditions for the transformation to be operationally 
unitary are therefore those of the special case discussed in section VI. The device has two computational input 
ports (labeled a and b) and four ancilla input ports (1-4). A projective measurement is made on four output 
ports (p,q,n,m) while the two remaining ports are the computational output (5 and 6). A photon of horizontal 
polarization represents a logical 0, and a vertically polarized photon represents a logical 1. We use the same 
notation as Pittman et al. [7]. For example, \H(V) a ) represents a horizontally(vertically) polarized photon 
in port 'a' and the Hadamard transformed modes are \F{S) a ) = \ [\H a ) ± \ V a )]. The four basis states of the 
computational input are |00) = \H a ) |if 6 _)_, |01) = \H a ) |V&)_, |10) = \V a ) |#b)_, |11) = \V a ) \V b ) and the output 
states are labeled as |00) = \H 5 ) \H 6 ) , |0l) = \H 5 ) \V 6 ) , |l0) = |V 5 ) \H 6 ) , |ll) = |K) \V 6 ) ■ The input ancilla 
state is 

IX) = \ (|#i) 1^4) \H 2 ) \H 3 ) + m \V 4 ) \H 2 ) \V 3 )) 
+ \ (\Vi) \H A ) \V 2 ) \V 3 ) + \V,) \V 4 ) \V 2 ) \H 3 )) , 

and the measurement projectors, P(l) = Iq® \ , represent the 16 possible success outcomes: 

\h) = \F p )\F g )\F n )\F m ) 

= \ (\H P ) + \V P )) (\H q ) + \V q )) (\H n ) + |K)) (\H m ) + \V m )) 
\Y 2 ) = \F p )\F q )\F n )\S m ) 

= \ (\H P ) + \V P }) (\H q ) + \V q )) (\H n ) + |K)) (\H m ) |K)) 



|*i5) - \S p )\S q )\S n )\F m ) 

= 1 (\H P ) |K)) (\H g ) \V q )) (\H n ) |K)) (\H m ) + |K)) 
\he) = \S P ) \S q ) \S n ) \S m ) 

= 1 (\H P ) |K)) (\H g ) \V q )) (\H n ) |K)) (\H m ) |K)) 

The polarizing beam splitters perform a unitary evolution on the input ports, characterized by the set of 
quantities U^. One can summarize the evolution of modes in He ® Ha to modes in Tic ® with the 
following linear map 



Hi) - 


- \H P ),\Vi)- 


>-i\V q ), 


H 2 ) - 


- \H 5 )AV 2 )- 


>\v 5 ), 


H 3 ) - 


- \H 6 ),\V 3 )- 


>\Ve), 


H A ) - 


+ \H m ),\V 4 )- 


- -i IK) 


H a ) - 


- \H q ),\V a )- 


» -i \V P ) , 


\H b ) - 


- \H n ),\V b )- 


>-i\V m ) 
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since ^* HlHp = 1, = — i, etc. As in the previous example, to evaluate the test operators, we first look at 

the terms 



\H q ) \H n ) 



u{\w)\x)) 

U{\Ql)\x)) 

u(\io)\x)) 
u{\n)\x)) 

The matrix elements of interest are now 



-i \H q ) \V m ) 



-i \V P ) \H n ) 
2 

- \V V ) \V m ) 



\H P ) \H m ) \H 5 ) \H 6 ) - i \H P ) \V n ) \H 5 ) \V 6 ) 
-i \V g ) \H m ) |Vb) |K) - \V q ) \ V„) |Vb) \H 6 ) 

' \H P ) \H m ) \H S ) \H 6 ) - i \H P ) |K> \H S ) |K) 
-i \V q ) \H m ) \V 5 ) \V 6 ) - |K) IK) |K) \H 6 ) 

\H P ) \H m ) \H 5 ) \H 6 ) - 1 \H p ) |K) \H 5 ) IK) ' 

-i |K) \H m ) IK) IK) - IK) \v n ) IH) \H 6 ) 

\H P ) \H m ) \H 5 ) \H 6 ) - 1 \H p ) |K) \H 6 ) |K) 

-i IK) \H m ) IK) IK) - IK) IK) IK) \H & ) 



((a\®(x\)Uip {L) U(\ X )l 
£ (M ® (xl) & \k£) \a) (a\ (k£\ U (\ X ) ® I/?)) 



(25) 



and the non-zero terms of the sum in (25) are 



\(m\(ki:\u(\x)^\oo))\ 2 = ^ 

|(0T|(fcZ|[/(| x )®|0T))| 2 = 1 

|(n|(fcZ|c/(| x )®|To))| 2 = 1 

|(To|(fcZ|[/(| x )®|n))| 2 = 1 



for all L. The test functions, {T(l)} , are then 



1 

64 



|ff„) |ff 6 > <JT 6 | (H a \ + \H a ) |K) (HI (H a \ 
+ |K) |^fe) (ff 6 | (K| + |K) |K) (HI (HI 



= m Isc 

and are indeed multiples of the unit operator in the computational input space. In this scheme T(£) = 1/64, and 
the probability of success is the sum of the individual probabilities of the 16 detection outcomes, Y^l t {l) = 1/4- 
The terms of the transformation matrices w^f can be calculated noting that the non-zero (k^ | (A| U \(3) \x) terms 



are 



(kE\ (oo|f/|oo)| x ) 

(K\ <0T|f/|01)| x ) 
(k£\ <TT| C/|10> |x> 
(kl\ <T0| C/|ll> |x> 



= e tipL - 3 



/8, 



where e 4< ^.° = l.e*^ 1 - 1 = -l.e^ 2 ' 1 = 1, 

a/3 



j e *<?i6,3 — 1 are phase factors of ±1. For this transformation to be 
operationally unitary, the vf^ matrices must all be proportional to each other. In certain outcomes, single-qubit 
operations (7r-phase shifts) are required to correct the phase factors so that the transformation is operationally 
unitary and the desired output is produced. The feed- forward processing matrices, Vfy, represent these single 
qubit operations. Setting 

T ->00,00 _ i0 LiO -,701,0! _ M, L1 t/TMT _ p i<Pl,2 t>T0,T0 _ w, ii3 
V (L) - 6 ' V (L) — C ' V (L) _ C ' V (L) — 6 ' 
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with all other elements equal to zero gives the appropriate corrections. The non-zero transformation matrix 
elements, are then 

00.00 , 

w L =1 

01,01 i 

w L =1 

11,10 i 

w L =1 
= 1 

for all L. Since the 16 evolution matrices are identical, the proportionality condition is satisfied. The transfor- 
mation is then 

0,1 

which is the CNOT operation. This gate takes the input state \tp) = ao 1 00) + ct\ 1 01) + a 2 1 10) + a s |11) and 
produces the state ao |00) + ai |0l) + a 2 |ll) + «3 |l0) with probability 1/4. Again, this map exhibits an 
effective nonlinear interaction between the photons since the formally equivalent unitary map is characterized 
by a nonlinear effective action operator Q (16). In this case one could choose 

Q = -y (3 + a|(l - h b ) + (1 - h b )a b )n a . 

Again, however, the operational unitarity is restricted to the subspace. Suppose we expand the computational 
subspace to include an extra photon in one of the input modes. As an example, consider the special state \S) = 
\H a ) \H b ) \H b ). The form of the projectors indicates that the detection events involve one and only one photon 
in the appropriate modes. Evaluating the corresponding test operator elements we find 

\(a\ (kZ\U{\ X )®\S))\ 2 = i), 

since the extra photon inhibits a success- indicating measurement result. The evolution cannot be operationally 
unitary in this expanded subspace because the test operator is no longer proportional to the unit operator. 



VIII. CONCLUSION 



In this paper we introduced a general approach to the investigation of conditional measurement devices. We 
considered an important class of optical iV-port devices, including those employing projectors of rank greater 
than unity, mixed input ancilla states, multiple success outcomes, and feed-forward processing. We also sketched 
how more general POVMs, rather than PVMs, could be included. The necessary and sufficient conditions for 
these devices to simulate unitary evolution have been derived. They are not surprising, and indeed from a 
physical point of view are fairly obvious. But to our knowledge they have not been discussed in this general 
way before. One of the conditions is that the probability of each successful outcome must be independent of 
the input density operator. Whether or not this holds can be checked by evaluating a set of test operators 
over the input computational Hilbert space, which is easily done for any proposed device. In the special case 
of only one successful outcome there is only one test operator to be computed; furthermore, if the ancilla state 
is pure and the success projector of rank one, then the passing of a test condition by that single test operator 
guarantees that the map is operationally unitary. In the case of more than one successful outcome it is a 
necessary consequence of operational unitarity that each of the test operators pass the test condition. This is 
not sufficient to imply operational unitarity in the multiple projector case unless the proportionality condition 
is also satisfied. The proportionality condition can often be satisfied by introducing feed-forward processing. 

Besides application in the analysis of particular proposed devices, we believe the general framework presented 
here will be useful in exploring the different types of pre-measurement evolution and measurements that might 
be useful in the design, optimization, and characterization of such devices. In particular, the conditional sign flip 
and polarization-encoded CNOT devices we considered functioned as operationally unitary maps only over the 
input computational subspaces for which they were originally proposed. So while effective photon nonlinearitics 
could be introduced, the degree to which they are physically meaningful is somewhat limited. An outstanding 
issue, perhaps even of interest more from the general perspective of nonlinear optics than from that of quantum 
computer design, is the study of potential devices that provide effective photon nonlinearities over much larger 
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input computational subspaces. The question remains: to what extent are such devices possible in theory and 
feasible in practice? 

Finally, we note that only in section VII did we assume that the pre-measurement unitary evolution U is 
associated with linear elements in an optical system. The more general framework of the earlier sections may find 
application in describing other proposed devices for quantum information processing that involve conditional 
measurement schemes in the presence of more complicated interactions [18]. 
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FIG. 1: A schematic diagram of a basic conditional measurement device. The input computational channels, He, and 
input ancilla channels, Ha, undergo unitary evolution. The measurement performed in the ancilla output space, Ha, 
indicates the success or failure of the computation. 
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FIG. 2: A schematic diagram of a conditional measurement device that incorporates feed-forward processing. The double 
line connecting the two small boxes represents a classical channel that carries the measurement result. Based on the 
outcome, the appropriate processing is performed on the output channel. 



